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Abstract 

Contributions of strange quarks to the mass and spin of the nucleon, charac- 
terized by the observables fj> s and As, respectively, are investigated within lattice 
QCD. The calculation employs a 2+1-flavor mixed-action lattice scheme, thus 
treating the strange quark degrees of freedom in dynamical fashion. Numerical 
results are obtained at three pion masses, = 495 MeV, 356 MeV, and 293 MeV, 
renormalized, and chirally extrapolated to the physical pion mass. The value ex- 
tracted for As at the physical pion mass in the MS scheme at a scale of 2 GeV 
is As = —0.031(17), whereas the strange quark contribution to the nucleon mass 
amounts to fx 3 = 0.046(11). In the employed mixed-action scheme, the nucleon 
valence quarks as well as the strange quarks entering the nucleon matrix elements 
which determine fx s and As are realized as domain wall fermions, propagators of 
which are evaluated in MILC 2+1-flavor dynamical asqtad quark ensembles. The 
use of domain wall fermions leads to mild renormalization behavior which proves 
especially advantageous in the extraction of fj< s . 

PACS: 12.38. Gc, 14.20. Dh 
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1 Introduction 

Strange quarks represent the lightest quark flavor not present in the valence component 
of the nucleon. Their study can thus provide insight into sea quark effects in the nu- 
cleon in isolated fashion. The two most fundamental properties of the nucleon are its 
mass and spin. The investigation presented here focuses on the extent to which those 
two properties are influenced by the strange quark degrees of freedom. The strange 
contributions to nucleon mass and spin can be characterized by the matrix elements 



respectively, where \N,j) denotes a nucleon state with spin polarized in the j-direction. 
In the case of the scalar matrix element, the vacuum expectation value, i.e., the vacuum 
strange scalar condensate, is subtracted; the intention is, of course, to measure the 
strangeness content of the nucleon relative to the vacuum. In the case of the axial matrix 
element, no subtraction is necessary since the corresponding vacuum expectation value 
vanishes. Note that As measures specifically the contribution of strange quark spin to 
nucleon spin; strange quark angular momentum constitutes a separate contribution not 
considered here. 

Aside from representing a fundamental characteristic of the nucleon in its own right, 
the scalar strange content /r s is also an important parameter in the context of dark 
matter searches [IrS]. Assuming that the coupling of dark matter to baryonic matter 
is mediated by the Higgs field, the dark matter detection rate depends sensitively on 
the quark scalar matrix elements in the nucleon, cf., e.g, the neutralino- nucleon scalar 
cross-section considered in [1] . One a priori reasonable scenario is that the strange quark 
furnishes a particularly favorable channel [Ij, since, on the one hand, it features a much 
larger Yukawa coupling to the Higgs field than the light quarks, and, on the other hand, 
is not so heavy as to be only negligibly represented in the nucleon's sea quark content. 
As a consequence, an accurate estimate of Jt s is instrumental in assessing the discovery 
potential for dark matter candidates. 

The contribution of strange quark spin to nucleon spin As is, in principle, more di- 
rectly accessible to experiment than fx s - As represents the first moment of the strange 
quark helicity distribution As(x) (including both quarks and antiquarks) as a function 
of the momentum fraction x. The helicity distribution can be determined via inclusive 
deep inelastic scattering and semi- inclusive deep inelastic scattering [5H7]. However, its 
extraction in practice still has to rely to a certain extent on assumptions about the 
dependence of As(x) on x, even in the semi-inclusive channels (which furnish direct 
information on As(x)), because of the limitations in accessing small x experimentally. 
Complementary information about As is obtained from the strange axial form factor of 




(1) 



and 




(2) 



2 



the nucleoli G S A (Q 2 ), which can be extracted by combining data from parity- violating 
elastic electron-proton scattering and elastic neutrino-proton scattering [8]. Extrapola- 
tion to zero momentum transfer, Q 2 = 0, again yields an estimate of As. Depending 
on the specific extrapolations and/or model assumptions adopted in determining As 
via the various aforementioned avenues, both significantly negative values for As have 
been put forward (SJE], as well as values compatible with zero [6J. An independent 
determination of As via lattice QCD, as undertaken in the present work, thus can be 
useful in several ways. Apart from shedding light on the fundamental question of the 
decomposition of nucleon spin, it can contribute constraints to phenomenological fits of 
polarized parton distribution functions. Furthermore, it influences spin- dependent dark 
matter cross sections [3J; although more accurate determinations of the scalar matrix 
elements discussed further above constitute the most urgent issue in reducing hadronic 
uncertainties in dark matter searches, As also plays a significant role in that context. 

A number of lattice QCD investigations of strange quark degrees of freedom in the 
nucleon have recently been undertaken [T0H26] . the majority of which have focused 
specifically on the scalar content. Studies of the latter have proceeded via two avenues: 
On the one hand, one can directly determine the matrix element (iV|ss|iV) via the 
appropriate disconnected three-point function; this methodology was adopted irfl [TQl 
[1 ^[T4T41T1 I19 [I22|, |2"4"] and also in the present work, as described in detail further below. A 
study of techniques suited to improve the efficiency of this approach has been presented 
in [27]. On the other hand, a somewhat less direct inference of the scalar strange quark 
content of the nucleon is possible via the study of the baryon spectrum, which is related 
via the Feynman-Hellmann theorem 

(BMB) = ^ (3) 

to the corresponding sigma terms for the baryon state \B) and quark flavor q. This 
avenue has been pursued in [T3 | [T5 | [2T | [23] . and a related methodology, combining lattice 
hadron spectrum data with chiral perturbation theory, was pursued in [25||26]. 

The characteristics of these various investigations of the scalar strange quark content 
of the nucleon are diverse. They include Nf = 2 calculations, in which the strange quark 
degrees of freedom are quenched [I2 HI3I23] . but also N f = 2+1 [TUlfToWnfTM2Tmil2ol2D] 
and even Nf = 2 + 1 + 1 [17,22j calculations. In some cases, lattice data at only one pion 
mass have been obtained to date and no extrapolation to the physical point has been 
attempted. The most stringent results obtained at the physical point including fully 
dynamical strange quarks were reported in [ToJ[T71|2B]. Ref. [T7J quotes m N f T Jm s = 
0.637(55) (74) in the N f = 2 + 1 case, and m N f T Jm a = 044(8) (5) in the N f = 2 + 1 + 1 
case; translated to fx a itself using m s = 95 MeV in the MS scheme at a scale of 2 GeV 
[28], these correspond to f Ta = 0.064(6) (7) for N f = 2 + 1 and f Ta = 0.044(8) (5) for 
Nf = 2 + 1 + 1. On the other hand, [15] and [26] report significantly lower values. In [15J, 



1 Strictly speaking, the method adopted in [T6l[T7] constitutes a hybrid of the two methods. 
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both direct three-point function calculations as well as indirect determinations via the 
Feynman-Hellmann theorem are used to arrive at a bound Jt 3 = 0.009(15)(16). In |26j, 
a more indirect analysis using lattice hadron spectrum data and chiral perturbation 
theory yields m^fTs = 21(6) MeV, which translates to fr s = 0.022(6). The results 
obtained in the present work are of comparable accuracy, fr a = 0.046(11), cf. ( 157)) . and 
are more consistent with the higher values reported in [T7] . 

The strange axial matrix element As, on the other hand, has also been investigated 
in [T0 ) [TT ) |24] . Apart from the exploratory study |10] , which, however, did not attempt to 
renormalize the results nor extrapolate them to the physical pion mass, these investiga- 
tions were based on dynamical quark ensembles containing only the two light flavors in 
the quark sea; the present lattice investigation, on the other hand, employs Nf = 2 + 1 
gauge ensembles, thus treating the strange quark degrees of freedom in dynamical fash- 
ion. The numerical results for As obtained on this basis are renormalized and chirally 
extrapolated, yielding the estimate As = —0.031(17), cf. (158)) . at the physical point in 
the MS scheme at a scale of 2 GeV. Within the uncertainties, this nevertheless remains 
compatible with the values obtained in the aforementioned other studies, though it is 
about 50% larger in magnitude. This suggests that systematic adjustments to the results 
quoted in those works through unquenching of the strange quark degrees of freedom, 
renormalization, and chiral extrapolation are not severe. 

Aside from the two quantities fx 3 and As considered in the present work, lattice 
investigations of the strange quark structure within the nucleon have also considered 
generalizations to non-zero momentum transfer, i.e., form factors [HI El], including 
calculations of the strange electric and magnetic form factors, which are of interest in 
the context of corresponding experimental efforts employing parity-violating electron- 
proton scattering [29]. Furthermore, also the strange quark momentum fraction and 
strange quark angular momenta in the nucleon have been investigated [20J . 

The present lattice investigation, a preliminary account of which was given in [SU] , is 
based on a mixed-action scheme developed and employed extensively by the LHP Collab- 
oration [511455] . The nucleon valence quarks as well as the strange quark fields appearing 
in the operator insertions in eqs. ([[]) and (j2J) are realized as domain wall fermions, propa- 
gators of which are evaluated in the background of (HYP-smeared) 2+1-flavor dynamical 
asqtad quark ensembles provided by the MILC Collaboration. Though computationally 
expensive, domain wall fermions lead to benign renormalization properties which prove 
especially advantageous in the case of the fo s observable; the substantial systematic un- 
certainties due to strong additive quark mass renormalizations encountered in analogous 
calculations using Wilson fermions [12l[2l] are avoided. The computational scheme is 
described in detail in section |2j Section [5J provides the raw numerical results, the renor- 
malization of which is discussed in section HI The renormalized results are extrapolated 
to the physical pion mass in section and systematic uncertainties and adjustments 
are considered in section (6J with conclusions presented in section [7J 
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2 Computational scheme 



2.1 Correlator ratios 

The lattice calculation of the nucleon matrix elements in (pQ), (JSJ) proceeds in standard 
fashion via correlator ratios of the type 



R[ r 



nuc -pofes 1 /_ 



^unpol 

a/3 



Eg Np{x,T)N a {0,0) 



(4) 



with nucleon interpolating fields N of the standard form (quoting here, for definiteness 
in flavor structure, the proton case) 



(5) 



where C denotes the charge conjugation matrix. The sums over spatial position x project 
the nucleon states onto zero momentum, whereas the sum over spatial position y is 
simply transcribed from (JXJ) , (J2J) . Since the nucleon contains no strange valence quarks, 
the three-point function averaged in the numerator of (jl]) factorizes, as written, into 
the nucleon two-point function and the strange quark loop. I.e., only the disconnected 
diagram, cf. Fig. [U contributes to the matrix elements under consideration. The strange 



N,t = T 



r,t = r 

O 







N, t = 



Figure 1: Disconnected contribution to nucleon matrix elements. The nucleon propa- 



gates between a source at t 
an intermediate time t = t. 



and a sink at t = T; the insertion of T 



^obs 



occurs at 



quark fields have already been reordered in the numerator of @ such as to make the 
standard minus sign associated with quark loops explicit. Finally, the T matrices allow 
one to choose the appropriate nucleon polarization and strange quark operator insertion 
structures. The denominator of @ corresponds to the unpolarized nucleon two-point 
function, obtained using 

1+74 



j-tunpol 



(6) 
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In the numerator, for the purpose of evaluating f Ts , unpolarized nucleon states are 
appropriate, corresponding to the choice T nuc = Y unpo1 ; furthermore, the scalar strange 
quark insertion is obtained by choosing T obs = 1, 

{R[T^ ol ,l}(r,T) - [VEV]) = R{f Ts } T ^T° f Ts (7) 



m N 

with the vacuum expectation value 

[VEV] = (-Z$ Sj (y,T)s 7 (y,T)) (8) 

to be subtracted. 

On the other hand, As is obtained by using the projector onto nucleon states polar- 
ized in the positive/negative j-direction in the numerator of (j4]), 

■pnttc -punpol 

as well as the operator insertion structure T obs = 77-75. Averaging over positive/negative 
j-direction (with a relative minus sign) as well as the three spatial j to improve statistics 
leads one to evaluate 

- i ■ 2 • \ Y>[ (-z 7i 75/2) T»** t 7 , 75 ](r, T) = R{As} T> ^° As , (10) 
6 3=1 

where the prefactor 2 compensates for the fact that the ratio (jlj) is normalized using the 
unpolarized nucleon two-point function in the denominator, even when the numerator is 
restricted to a particular polarization. Lastly, the prefactor (— i) cancels the additional 
factor % which the 75-matrix in the operator insertion acquires when the calculation is 
cast in terms of Euclidean lattice correlators; thus, the As obtained through ffTUj) is 
already Wick-rotated back to Minkowski space-time. Note that ffTUl) does not call for 
the subtraction of a vacuum expectation value, since the latter vanishes in the ensemble 
average. In practice, this numerical zero was nevertheless subtracted from (TlOT) in order 
to reduce statistical fluctuations. 



2.2 Lattice setup 

The averages in the correlator ratios (jlj) were carried out using the three MILC 2 + 1- 
flavor dynamical asqtad quark ensembles listed in Table [TJ HYP-smearing was applied to 
the configurations. Both the valence quarks in the nucleon two-point functions and the 
strange quark fields appearing in the matrix elements (jl]) and (j2j) were implemented using 
the domain wall fermion discretization, with parameters L5 = 16, M5 = 1.7. The domain 
wall quark masses, also listed in Table (U are fixed by the requirement of reproducing 
the pion masses corresponding to the MILC ensembles [31]. This mixed action setup has 
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asq 

am l 


am a s sq 


am}> w *' 




m n 


# configs 


m N 


0.007 


0.05 


0.0081 


0.081 


293 MeV 


468 


1107 MeV 


0.01 


0.05 


0.0138 


0.081 


356 MeV 


448 


1155 MeV 


0.02 


0.05 


0.0313 


0.081 


495 MeV 


486 


1288 MeV 



Table 1: N f = 2 + 1, 20 3 x 64 MILC asqtad ensembles with lattice spacing a = 0.124fm 
used in the present investigation. Uncertainties in the pion and nucleon masses extracted 
from the corresponding two-point functions are under 1%. 

been employed extensively for studies of hadron structure by LHPC [31j[33] and further 
details concerning its tuning can be found in the mentioned references. 

The space-time layout of the calculation is shown in Fig. |2J The strange quark loop 
trace was evaluated using stochastic estimation. Positioning the nucleon source at time 
t = 0, bulk complex Z(2) stochastic sources with support in all of space within the tem- 
poral range t — 3, ... ,7 were introduced. This corresponds to averaging the correlator 
ratio (jl]) with respect to the operator insertion time r over the aforementioned time 
slices (after having duly divided out the length of the temporal range). The stochastic 
estimate of the strange quark loop trace was performed employing 1200 of the described 
stochastic source vectors per gauge sample. In particular, obtaining a signal for As 
depends on accumulating high statistics in the stochastic estimator. The scalar matrix 
element, on the other hand, requires less statistics in the strange quark loop trace, but is 
more susceptible to gauge fluctuations. The sink time T at which the nucleon two-point 
functions are contracted and projected onto zero momentum remains variable in this 
layout. 

While this scheme provides a high amount of averaging in relation to the number 
of strange quark propagator inversions, it precludes testing for a plateau in the three- 
point function by varying the operator insertion time r. The positioning of the bulk 
stochastic source was motivated by previous work [3TJ[33] in the same mixed-action 
scheme, which investigated (connected contributions to) a wide variety of observables, 
and the results of which indicate that, having evolved in time from t = to t = 3, 
excited state contaminations in the system are already small compared to the statistical 
uncertainty of the calculation performed here. To achieve this suppression of excited 
state contributions, the quark fields in the nucleon sources fl5]) are Wuppertal-smeared 
such as to optimize the overlap with the nucleon ground state [3"TH3"3"] . Despite the 
relative positioning of the nucleon source and the strange quark insertion time range 
being fixed in this manner, it does nevertheless remain possible to exhibit the extent 
to which results depend on the separation between operator insertion and nucleon sink 
time T, which is still variable in the present scheme. Correlator ratios will be shown 
further below as a function of T, with asymptotic behavior being seen to emerge three 
lattice time steps beyond the end of the operator insertion range, i.e., for sink times 
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— I 1 1 > t 

3 7 

Figure 2: Setup of the lattice calculation. The nucleon source is located at lattice time 
t = 0. An average over operator insertion times r is performed in the range r = 3, . . . , 7; 
accordingly, stochastic sources are distributed over the bulk of the lattice in this entire 
time range. The temporal position T of the nucleon sink is variable. 

T > 10. This corroborates the suitability of the choice for the bulk stochastic source 
time range adopted in the present calculation. 

The correlator ratio (J4]) exhibits statistical fluctuations not only due to the strange 
quark loop factor discussed above, but also due to the nucleon two-point function factors. 
To reduce these fluctuations, it is equally necessary to sample the latter to a sufficient 
extent. To this end, multiple samples of the nucleon two-point function were obtained 
by employing eight different spatial positions for the nucleon source on the source time 
slice. In addition to this eight-fold sampling of the nucleon two-point function, the de- 
scribed scheme can be accommodated several times in temporally well-separated regions 
on the lattice; in practice, three replicas of the entire setup specified above were placed 
on the lattice, separated by 16 lattice spacings in the time direction, thus further en- 
hancing statistics. Each gauge configuration therefore yielded altogether 24 samples of 
the numerator and denominator in the correlator ratio Note again that averaging 
was further improved in the case of As by taking into account nucleon polarization axes 
aligned with all three coordinate already made explicit in ( TTOj) . 

3 Numerical results 

The numerical results for the correlator ratios R{fr a } an d i?{As}, cf. (JTj) and (fT0~|) . 
averaged over the insertion time r as described in section I2.2[ are shown in Figs. |3] and 
S] as a function of sink time T. The correlator ratios start out near vanishing values at 
small sink times T, and then accumulate strength as T traverses the region of support 
of the stochastic strange quark source. Beyond this region, R{fT s } and i?{As} begin to 
level off and approach their asymptotic values Jt s and As, respectively. There is only 
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a limited time window in which the latter behavior can be observed due to the increase 
in statistical fluctuations for large T. 





R{As}\ T=w 


fl{As}| T =10,...,12 


293 MeV 


-0.023(17) 


-0.023(25) 


356 MeV 


-0.030(9) 


-0.030(12) 


495 MeV 


-0.021(4) 


-0.022(5) 



Table 2: Correlator ratio R{As} at sink time T = 10, and averaged over T = 10, . . . , 12. 
In the latter case, the error estimate is obtained by the jackknife method, i.e., the 
correlations between values of i?{As} at neighboring sink times are taken into account. 

The correlator ratio R{As} consistently behaves in line with the expectation ad- 
vanced in section 12.21 namely, there is no significant trend in the correlator ratio beyond 
T = 10. This is quantified in Table EJ which compares the value of _R{As} at T = 10 
with its T = 10, . . . , 12 average. 

The behavior of the lattice data for the correlator ratio R{fr s } is not as smooth as in 
the case of i?{As}. In particular, in the m f = 356 MeV correlator ratio, one notices an 
enhancement in the T = 12, . . . , 14 region compared with the value at T = 10. On the 
other hand, both in the m n = 495 MeV and the = 293 MeV data, plateaux appear 
to be reached at T = 10. Table [3] again provides a comparison of the T = 10 value of 
R{fr s } with its T = 10, . . . , 12 average; in the case of = 356 MeV, the results are 
barely compatible within the statistical uncertainties. Nevertheless, an interpretation 
of this behavior as a statistical fluctuation, as opposed to a systematic effect, seems 
most plausible: The direction of the deviations from the expected plateau behavior 
in the correlator ratio beyond T = 10 is not consistent across the data sets, with 
R{fr s } decreasing slightly for m n = 495 MeV, while it rises in the other two cases. 
Furthermore, there is no clear trend as a function of pion mass, with the largest upward 
deviation occurring for the middle pion mass, m n = 356 MeV, while at = 293 MeV 
the deviation is again quite small. 

In summary, thus, at the level of statistical uncertainty achieved in the present 
calculation, no systematic excited state effects stemming from a too restricted source- 
sink separation can be reliably extracted; or, in other words, such effects are smaller 
than the aforementioned statistical uncertainties. In view of this, in the following, the 
T = 10 values of the correlator ratios i?{As} and R{fr a }, as reported in Tables [2] and EJ 
will be regarded as the most reliable estimates for their asymptotic limits As and fx s , 
respectively. Systematic uncertainties, including the ones due to excited states, will be 
revisited and discussed in detail in section [6j 
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gure 3: Correlator ratio R{fr a }, cf. flZ}, averaged over insertion time r as described 
section \'2.'2\ as a function of sink time T, for the three pion masses considered. 
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Figure 4: Correlator ratio R{As}, cf. ffTUl) . averaged over insertion time r as described 
in section I2.21 as a function of sink time T, for the three pion masses considered. 
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R{fT 3 }\T=lO 


-^{/t s }|t=10,...,12 


293 MeV 


0.060(10) 


0.064(14) 


356 MeV 


0.043(6) 


0.051(8) 


495 MeV 


0.046(3) 


0.045(4) 



Table 3: Correlator ratio R{fr s } at sink time T = 10, and averaged over T = 10, . . . , 12. 
In the latter case, the error estimate is obtained by the jackknife method, i.e., the 
correlations between values of R{fr s } at neighboring sink times are taken into account. 

4 Renormalization 

To establish a connection with phenomenology, quantities measured on the lattice need 
to be related to their counterparts in a standard renormalization scheme such as the MS 
scheme at a scale of 2 GeV. An advantage of the domain wall fermion discretization, used 
in the present treatment to represent the strange quark fields entering the quark bilinear 
operator insertions in the matrix elements (JT|) and (T5]), lies in its good chiral symmetry 
properties, which lead to a mild renormalization behavior. As a result, even though 
not all elements necessary for a complete renormalization of the quantities considered 
are available (e.g., matrix elements of gluonic operators in the nucleon, with which the 
flavor singlet parts of the strange quark bilinears may mix), it is possible to estimate the 
associated uncertainties in the renormalization, indicating that these do not dominate 
over the statistical uncertainties of the calculation. 

4.1 Scalar matrix element 

Consider first the renormalization of Jt s , specifically, the proposition that the combi- 
nation m s (N\ss\N) can be treated as independent of the renormalization scheme and 
scale, 

m s (N\ss\N)\ ienoim = m s (N\ss\N)\ hare . (11) 

As has been emphasized previously in [14J, this behavior is contingent upon chiral sym- 
metry being maintained in the lattice evaluation of the bare quantities. In general, 
the strange scalar operator ss can mix both with light quark and gluonic operators; 
decomposing it into flavor singlet and octet components, 

(12) 

(where q G {u, d, s}), the flavor singlet part can mix with the relevant gluonic operator, 

(qq) Teaoim = Zf{qq) h&ve + Zf(F 2 ) h&ve (13) 



l 

* S= 3 



qq 



V3q\t 
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and also the renormalization constant Z| 8 of the octet part, 

(gA 8 g) r ™ = Zf(q\ s q) h ™ e (14) 

in general is not identical to Z®P. However, these mixing effects would have to proceed 
[H] via diagrams in which the qq insertion occurs within a disconnected quark loop; 
furthermore, since qq = qiqR + qnqL, a chirality flip would have to take place within 
that loop. Gluon vertices coupling to the loop cannot effect such a chirality flip. Thus, 
if all explicit chiral symmetry breaking effects are excluded, both in the dynamics as 
well as by adopting a mass-independent renormalization scheme, the gluonic admixture 
in ffTB"]) is avoided. By extension [T4] , also no mixing of ss with uu + dd can take place, 
i.e., the renormalization factors of the flavor singlet and octet parts in ( !T2|) are equal 
under the assumption of strict chiral symmetry. Then, ss renormalizes in a purely 
multiplicative manner. 

Complementarily, also the strange quark mass m s only renormalizes multiplicatively, 
m renorm _ z m m^ are , when chiral symmetry is maintained. For generic non-chiral lattice 
fermion discretizations, the renormalized strange quark mass 777^ enorm acquires strong ad- 
ditive contributions supplementing the aforementioned multiplicative renormalization, 
invalidating the behavior (fTTj) . This is again avoided when chiral symmetry is respected. 
Indeed, when renormalization is purely multiplicative, in view of the Feynman-Hellmann 
theorem, 

m s (N\ss\N) = m s ^ , (15) 

factors Z m acquired by the strange quark mass under renormalization cancel on the 
right-hand side of f lT5|) . implying that also the left-hand side is invariant. 

Chiral symmetry of the lattice fermion discretization is therefore instrumental in 
establishing a simple renormalization behavior of fx a , and, indeed, is realized to a good 
approximation in the present treatment due to the use of domain wall fermions in the 
coupling to the operator ss. Whereas exact chiral symmetry would only be achieved 
using an infinite fifth dimension in the domain wall construction, in practice, the extent 
of this dimension has been chosen sufficiently large as to render the residual mass m res 
quantifying the violation of chiral symmetry an order of magnitude smaller than the 
light quark mass in all ensembles considered [31]. One must, however, be careful in 
assessing the significance of this smallness, since, on the other hand, the light quark and 
gluonic operators with which ss mixes may have nucleon matrix elements much larger 
than (iV|ss|./V); for (N\uu + dd\N) , this is certainly the case. This can potentially offset 
the suppression of m res . One may speculate that mixing with the operator uu + dd 
constitutes the strongest effect, since the light quark fields are special in that they 
form the valence component of the nucleon, which has no counterpart in the vacuum 
expectation value that is subtracted off throughout, cf. (pQ). On the other hand, the 
presence of the valence quarks also strongly distorts the gluon field in the nucleon. No 
estimate of the gluonic admixture to ss is available, but the light quark admixture under 
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renormalization will be argued below to constitute an effect of the order of 1%. In view 
of the statistical uncertainties associated with the present determination of /y s , which 
amount to about 20%, a putative gluonic mixing effect would have to be an order of 
magnitude larger than the light quark mixing effect in order to appreciably influence 
the final result for fx s - This seems a rather implausible scenario. For this reason, the 
strong constraint on mixing with light quarks derived below will be taken as indication 
that violations of f fTTj) are negligible at the present level of statistical accuracy of fc a ■ 

Concentrating thus on the effect stemming from the mixing with the operator uu+dd, 
an estimate of the possible residual violation of (llip can be obtained from the Feynman- 
Hellmann theorem as follows [14|. The residual breaking of chiral symmetry can be 
parametrized to leading order via the additive mass renormalization m reS (J , which in 
general depends on the flavor q for which one is considering the domain wall Dirac 
operator, 

m~ = Z m (m h q ^ + m res , q ) . (16) 

Using ( TT5|) . one has 



, renorm 



m S (N\ss\N)\ haxe = ^E ^Zm^oare ("0 



q 



« m s (N\ss\N)\ lenoim + m~^l(N\uu + dd\N)~ 

where m res j denotes the light quark flavor residual mass, and where only the dominant 
correction has been kept in the second line by using m^ are 3> m reStS and (N\uu+dd\N) 3> 
(N\ss\N). No direct data for the factor dm res j/dm^ a,Te are available, but a rough order of 
magnitude estimate can be constructed from related data on the residual mass obtained 
within the LHPC program and in the present work. This estimate is discussed in the 
Appendix, with the result 



dm. 



resl 



w -0.00035 

am a =0.081 , am; =0.0081 



cf. (j42|) . at the lowest light quark mass considered in the numerical calculations in this 
work. In turn, an upper bound for the light quark scalar matrix element (N\uu + 
<i<i|iV) renorm is given by its magnitude at the physical point; combining a typical phe- 
nomenological value for the nucleon sigma term, mi{N\uu + dd\N) pa 60MeV, cf., 
e.g., [31] for a recent discussion, with the physical light quark mass |28j, mi 3.5 MeV, 
yields (N\uu + dd\N) ~ 17. Supplementing this with the value of the strange quark 
mass [28] , m s = 95 MeV (all aforementioned quark masses being quoted in the MS 
scheme at 2 GeV) , yields 

5(m s (N\ss\N)) Pa -0.6MeV (19) 

for the correction term in the second line of ( TITI) . This would amount to a 1% upward 
shift of the bare result for /y s as one translates the quantity to the MS scheme at 2 GeV. 
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Compared to the statistical uncertainties of the present calculation, a correction of this 
magnitude is negligible, and in the following, ( I lip will therefore be taken to hold for the 
present calculation within its uncertainties^]. 

It should again be noted that the mixing with gluonic operators has not been quan- 
tified in the above considerations. Whereas the weakness of the mixing with light quark 
operators (which, after all, is mediated by the coupling to the gluonic fields) makes it 
seem improbable that mixing with gluonic operators themselves is significant compared 
to other uncertainties of the calculation, explicit corroboration of this expectation would 
be desirable. 



4.2 Axial vector matrix element 

Turning to the axial vector matrix element, chiral symmetry again provides important 
constraints, although it cannot completely exclude mixing effects, due to its anoma- 
lous Ua(1) breaking. The domain wall fermion discretization admits a five-dimensional 
partially conserved axial vector current A ^ obeying a Ward-Takahashi identity of the 
form [SSJEB] 

A,Al = 2mJ« + 2J« g (20) 

The first term on the right-hand side represents the explicit breaking of chiral symmetry 
by the quark masses, whereas the second term in the flavor-octet case encodes the 
residual chiral symmetry breaking present for a fifth dimension of finite extent [36J, 
Jg w m Tes . Thus, up to these chiral symmetry breaking effects, which will be revisited 
further below, the flavor octet part of the current A't is conserved and undergoes no 
renormalization. On the other hand, in the flavor singlet case, the second term on 
the right-hand side of ( 120]) in addition encodes the coupling to the gluonic topological 
density which leads to the anomalous Ua(1) breaking of chiral symmetry. This opens 
the possibility of operator mixing under renormalization in the flavor singlet component. 
With respect to direct gluonic admixtures, the axial case differs somewhat from the 
scalar case discussed further above. The continuum axial vector current A^° nt ^ receives 
no direct gluonic admixtures [37], since the relevant gluonic operator would be the 
operator which, upon taking the divergence, yields the gluonic topological density, 
dfxK^ = g 2 FF. However, itself is not gauge invariant, and therefore the gauge 
invariant operator A^ ont ^ receives no admixtures from K^. While it is not clear to what 
extent this argument is modified at finite lattice spacing, the fact that direct gluonic 
admixtures to the flavor- singlet axial vector current must vanish in the continuum limit 
suggests that any such modifications would be small compared to the mixing of the 
strange quark axial current with the light quark axial currents. Concentrating thus on 

2 Note that the consideration of the factor dm res ^i/ drn)^ 11 ™ given in the Appendix relies on a number 
of assumptions and should be viewed as no more than a rough estimate of what is most likely an upper 
bound on the magnitude of this correction factor. In view of this, applying the correction (|19|) to the 
numerical data does not seem warranted. 
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the latter, one can obtain a corresponding estimate of operator mixing effects as follows. 
Decomposing 

= g qiulM ~ ^3g7/,7sA 8 g (21) 

(where q G {u, d, s}), and allowing for different renormalization constants for the singlet 
and octet parts, 

/ — \ renorm ^ ^00 — 

(S7 M 75S) = -Z^q^q - —Z^q^ 5 X 8 q (22) 

1 700 _ 788 

= Zfs W s + L - Za Za Zfq^q (23) 

The relative strength of mixing — Z^)/Z^ is not available for the specific lat- 
tice scheme used in the present work, but has been estimated for the case of clover 
fermions [TTIEB]. To ^ ne ex t en t that this encodes the effect of the axial anomaly as 
opposed to specific lattice discretization effects, it can be taken as indicative of the 
strength of mixing also in other lattice schemes such as the one used in the present in- 
vestigation. Inasfar as it is influenced by the lattice scheme, it presumably can be taken 
to represent an upper bound for the strength of mixing in schemes which better respect 
chiral symmetry such as the one used here (assuming there are no accidental cancella- 
tions). Quantitatively, for the conversion into the MS scheme at the scale /i = 2 GeV, 
one obtains in the two-flavor clover fermion case [11] 

Zf-Zf 0.0082 , . 

A gg A = = 0.011 24 

Zf 0.765 V ; 

Correcting by a factor 3/2 to translate to the three-flavor case considered here, and 
supplementing with Z^ 8 1.1, as obtained below, as well as A(m + d) 0.42 (at the 
physical pion mass, cf. [33]), the uncertainty from mixing with light quarks, i.e., the 
second term in fl23l) is estimated to amount to 



S(As) w 0.0025 (25) 

directly at the physical pion mass. Due to the indirect nature of this estimate, the shift 
(|25|) will not be applied to the (chirally extrapolated) central value for As obtained in 
this work, but will be treated as a systematic uncertainty. 

It should be remarked that the considerations for the clover case [UJ referred to 
above pertain to lattices with a substantially finer spacing than used here (0.073 fm 
vs. 0.124fm). However, adjusting for this is expected to modify (1251) merely by a few 
percent and thus is negligible in the present context. This can be inferred from the 
magnitude of the 0(a) improvement corrections to the renormalization constants quoted 
in [UJ; note also that the present HYP-smeared mixed action scheme, which is fully 
0(a)-improved, suffers only from very benign finite lattice spacing effects, as evidenced, 
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e.g., by the congruence between nucleon mass measurements in this same scheme and 
corresponding MILC determinations on much finer, a = 0.06 fm lattices [32J. 

Turning to the renormalization constant Z^ 8 , the axial vector quark bilinear used 
in practice in evaluating As is the local s(a;)7 M 753(2) as opposed to the corresponding 
flavor component of the partially conserved of eq. ( 120]) . To renormalize this local 
operator, the standard scheme [SB] can be applied, with the modification that, in the 
present case, it is not the pion current but the 77 current which is relevant: Calculating 
the (connected parts of the) two-point functions of both the conserved and the local 
current, 

C(i + 1/2) = ^2(J§(x,t) 77(0,0)) | conn 

X 

L{t) = ^(g 7 o75(A 8 /2)gr/(0,0))| conn , (26) 

X 

Z^/Z^ can be extracted from an appropriate ratio which takes into account the tem- 
poral offset between the two currents, 

1 ( C(t + 1/2) + C(t- 1/2) 2C(t + l/2) \ t/a>i Zf_ 

2 V 2L(t) L{t) + Lit + 1)) Z™ 1 ' 

Note that the full correlators (?7o7s(A8/2)g 77(0, 0)) and (Aq(x, t) 77(0,0)) acquire dis- 
connected contributions for unequal strange and light quark masses; however, for the 
specific purpose of extracting Z^/Z^, any ratio of quantities which only differ by this 
overall renormalization factor is suitable, including using only the connected parts of 
the aforementioned correlators, as indicated in (j2"6"|) . Table H] lists the values obtained 
for Z^/Z^, which are applied to the bare lattice measurements of As. 



m bare 


0.0081 


0.0138 


0.0313 


0.081 


788 / 788 
I^A 


1.09 


1.09 


1.10 


1.13 



Table 4: Renormalization factor Z^/Z^ at varying m} 5are . 

A systematic uncertainty is associated with this scheme of renormalizing As, due to 
residual sources of chiral symmetry breaking. One of these sources is the finite extent 
of the fifth dimension in the domain wall fermion construction. At finite m res , the 
renormalization constant of the partially conserved current, Z^ 8 , can deviate from the 
unit value it would take if chiral symmetry were strictly observed [39]. Secondly, note 
that, at finite lattice spacing a, there is a certain tension between adopting a mass- 
independent lattice renormalization scheme and maintaining 0(a) improvement [40J. 
The renormalization constant Z^ in general contains dependences of order 0(m q a), 
evident in the slight variation displayed in Table HJ Since the lattice data necessary to 
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perform the continuum limit a — > are not available, two options remain: One option 
would be to extrapolate Z^/Z^ to the chiral limit, thus obtaining a mass-independent 
renormalization scheme in more direct correspondence to the MS scheme, but spoiling 
0(a) improvement. On the other hand, by retaining the leading quark mass dependence, 
i.e., applying the finite-m g renormalization constants in Table S] ensemble by ensemble, 
0(a) improvement is maintained, at the expense of introducing a slight mass dependence 
into the renormalization scheme at finite a. The mass dependence, implying a breaking 
of chiral symmetry in addition to the one encoded in the residual mass m res , is then 
expected to be of order 0(m q a 2 ). Since the present investigation yielded lattice data at 
only a single lattice spacing a, precluding a direct estimate of the effects of finite lattice 
spacing, maintaining 0(a) improvement seems sufficiently desirable to elect the latter 
alternative, i.e., applying the renormalization constants in Table H] ensemble by ensemble 
and thus introducing a slight mass dependence into the renormalization scheme. 

A way of estimating the systematic uncertainty in the renormalization of As re- 
sulting from the residual breaking of chiral symmetry due to the above sources lies in 
the mismatch between the axial vector and the vector renormalization factors, 
vs. Zy/Zy, which would remain equal if chiral symmetry were strictly maintained |¥Tj . 
In the lattice scheme used here, these factors typically differ by 3% [33JE An additional 
systematic uncertainty of this magnitude will therefore be attached to the renormalized 
value of As. 



5 Chiral extrapolation 

Chiral extrapolation formulae for strange quark matrix elements in the nucleon have 
been given in [32] • At leading order (LO), both (JV|ss|iV) and As are constant in m\. 
On the other hand, when one evaluates one-loop effects, including A-resonance degrees 
of freedom, one obtains next-to-next-to-leading-order (NNLO) formulae which contain 
too many parameters to be effectively constrained by the restricted set of lattice data 
at three pion masses accessed in this work. However, practicable fits can be constructed 
by reducing the NNLO formulae given in [32] to the chiral effective theory without A- 
resonance degrees of freedom, which is achieved by setting = and correspondingly 
also eliminating the counterterms associated with the A-resonance degrees of freedom. 
In this case, the behavior of (iV|ss|iV) reduces to a linear function in the light quark 
mass, i.e., in m^, 

(N\ss\N) = Si + S 2 ml (28) 
with the two fit parameters S\ and 5*2, whereas As retains a chiral logarithm, 



As = D l 



^9a 2 i / 2 / 2\ 



+ D 2 mi (29) 



3 To be precise, [33J considered the isovector currents, i.e., determined /Z^ and Z^p /Z^ 3 . 
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with the two fit parameters D x and D 2 ; the dependence on the scale // is of course 
absorbed by the D 2 counterterm. The pion decay constant / is normalized such that 
/ ~ 132 MeV, and the physical axial coupling constant is qa ~ 1.26. Corresponding 
chiral fits to the renormalized lattice data for m s (N\ss\N) = m^fTs and As are shown 
in Fig. [5j The LO constant fits and the (reduced) NNLO fits are consistent with one 
another. Positing the LO constant behavior in leads to artificially low estimates of 
the uncertainties; in this case, the statistical error bars are dominated by the most ac- 
curately determined = 495 MeV data points, at which, on the other hand, the chiral 
extrapolation formulae are the least trustworthy. Plausible estimates of the statistical 
uncertainties of the extrapolated values are given by the (reduced) NNLO fits, which 
allow for variation of the observables with m^. The estimates of fx a and As resulting 
from the (reduced) NNLO fits at the physical pion mass in the MS scheme at a scale 
of 2 GeV are (before taking into account systematic effects, which are discussed in the 
next section), 

f Tt = 0.057(11) , (30) 

where the physical nucleon mass has been used to convert the fit result from Fig. [5] back 
to f Tm , and 

As = -0.028(14) . (31) 

Note again that the quoted uncertainties at this point contain only the statistical error 
from the lattice measurement, propagated through the chiral extrapolation; systematic 
uncertainties and adjustments of the results (I3~0~I) and (I3II) are elaborated upon in the 
next section. 

6 Systematic corrections and uncertainties 

Several sources of systematic uncertainty should be taken into account with respect to 
the two results ( 1301) and (fSTT) : 

Renormalization uncertainties: Uncertainties associated with renormalization were 
already discussed in section HI In the case of fx a , uncertainties due to mixing with light 
quark operators generated by residual breaking of chiral symmetry were estimated to 
be of the order of 1% and will thus not be considered further here. In the case of As, 
mixing effects were less well constrained because of the anomalous breaking of chiral 
symmetry. The potential correction to As was estimated to amount to 5(As) ~ 0.0025 
towards less negative values, cf. (l2"5]) . In addition, a 3% uncertainty was assigned to the 
renormalization factor Z A 8 due to residual sources of chiral symmetry breaking. These 
included finite-m res effects as well as effects of adopting a not fully mass-independent 
lattice renormalization scheme, in order to preserve 0(a) improvement. Adding the 
uncertainty from operator mixing and the one associated with Z A 8 in quadrature implies 
an uncertainty in As of +0.003/ — 0.001. 
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Figure 5: Pion mass dependence of the results for m s (N\ss\N) = m^f^ an d As. 
Filled circles represent renormalized lattice data; in the case of As, these are obtained 
by multiplying the T = 10 values from Table |2] by the corresponding renormalization 
constants from Table HI whereas mjv/r s is obtained by multiplying the T = 10 values 
from Table [3] by the corresponding nucleon masses from Table [lj Open symbols show 
chiral extrapolations of the lattice data to the physical pion mass, cf. main text. Open 
circles represent the LO (constant) chiral extrapolations, whereas open squares represent 
the reduced NNLO extrapolations obtained by dropping the A-resonance degrees of 
freedom, with the dashed lines showing the pion mass dependences of the central values 
in the latter case. 
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Finite lattice spacing effects: Since the present investigation only employed ensem- 
bles at a single lattice spacing, a = 0.124 fm, no direct assessment of the a-dependence 
of the results was possible. This motivated the insistence on a fully 0(a)-improved cal- 
culational scheme, in order to minimize the influence of the finite lattice spacing from 
the outset (and, in the process, part of the lattice spacing dependence was already sub- 
sumed under the uncertainty in renormalization, as noted above). In the case of Jt s , 
which is related via the Feynman-Hellmann theorem to the nucleon mass, an estimate 
of the uncertainty due to discretization effects can be inferred from the a-dependence 
of the nucleon mass. As shown in [32], already at a = 0.124fm, the nucleon mass in 
the present HYP-smeared mixed action scheme coincides with the MILC a = 0.06 fm 
results. Given that the MILC results themselves still change by about 10% going from 
a = 0.124fm to a = 0.06 fm, the residual 0(a 2 ) effect in the latter case is expected to 
be around 3%. This will therefore be taken as the generic estimate of the magnitude of 
finite lattice spacing effects in the present calculation, both for fa and As. 

Uncertainty due to truncation of the chiral perturbation series: Taking the 
deviation between the LO and the (reduced) NNLO fits at the physical pion mass, 
cf. Fig. [5j as a measure of the uncertainty due to truncation of the chiral perturbation 
series, a 6% uncertainty is attached to the value of Jt s and a 14% uncertainty to As. 

Effects of inadequate lattice dimensions: Both the spatial extent of the lattice 
and the temporal separations between nucleon source, sink and operator insertion are 
limited. Consequently, results are influenced both by interactions with periodic copies of 
the lattice as well as excited state admixtures. Neither of the corresponding uncertainties 
were directly quantifiable within the present calculation. On the one hand, ensembles 
with only a single lattice extent were employed; on the other hand, no systematic excited 
state effects could be gleaned from the sink position dependence of the lattice data at 
the present level of statistical accuracy, as discussed in section [3j An indication of 
the possible magnitude of such effects can be inferred from lattice calculations of the 
nucleon axial charge g^, the isovector light quark analogue of As, which represents a well- 
studied benchmark quantity. Calculations of g& within the present lattice scheme [33] 
and others [USUI] exhibit a deviation from the phenomenological expectation of up to 
10%, with the cause of this deviation attributed to either excited state contaminations 
or finite lattice size effects. This will therefore be taken as an estimate of the systematic 
uncertainty due to such effects also in the present calculation. 

Adjustment of the strange quark mass: The strange quark mass am a s sq = 0.05 in 
the gauge ensembles used in the present calculation lies appreciably above the physical 
strange quark mass, which a posteriori was determined to be am a s sq ' phys = 0.036 [16J. 
For the case of the strange scalar matrix element, a corresponding correction factor was 
estimated in [16], namely, 

jA^^% = -2-2 ■ 0.31 fm = -0.68 fm (32) 
om s orris 
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Multiplying this by the shift in m a s sq , 8m a s sq = (0.036 — 0.05)/a, yields, in view of the 
Feynman-Hellmann theorem, an enhancement of the strange scalar matrix element by 

6((N\ss\N)r« = 6m?« ®^ a N sg = 0.077 (33) 

om s orris 

To translate this to the present scheme, one needs to rescale the strange quark mass, 
m asq _ (jff-sq /rn DWF )mf . The rescaling factor m asq /m DWF varies only weakly be- 
tween am a s sq = 0.05 and am a s sq = 0.02, namely, from 0.617 in the former case to 0.639 
in the latter, cf. Table [TJ Interpolating linearly, the most appropriate value for the 
present consideration is the one halfway between am a s sq = 0.05 and am a s sq = 0.036, i.e., 
m asq /m DWF = 0.622. Since (155|) has two powers of m a s sq in the denominator and one in 
the numerator, altogether, translated to the present scheme, 

5((N\ss\N)) DWF = -^SiiNMN)) " = 0-048 (34) 

To obtain a measure of the relative change in (iV|ss|iV) implied by this, note that the 
result d30D, multiplied by m N /m° WF = 7.3, yields (N\ss\N)^j^ F= ^ mQ81 = 0.42, and 
therefore the adjustment (1341 corresponds to an enhancement of (iV|ss|iV) by a factor 
1.115 as one shifts the strange quark mass to its physical value. 

On the other hand, Jt s itself acquires an additional factor m DWF,phys j m DWF ^ 
m asq,phys J m asq _ q j2 as one shifts the strange quark mass to the physical point, im- 
plying that fx s is reduced by a factor 0.72 • 1.115 = 0.80. The reduction of m s in fact 
overcompensates the enhancement of (JV|ss|iV). Altogether, thus, (150]) will be adjusted 
by a factor 0.80 to arrive at the physical value; in addition, a systematic uncertainty of 
3% will be associated with that adjustment in view of a 15% uncertainty in ( 1521 . cf. [16J, 
as well as the variability in m asq /m DWF . 

For the case of As, no similarly detailed consideration is available. However, it seems 
plausible that the leading effect of a lowering of the strange quark mass is an overall 
enhancement of the strange quark density, with the detailed dynamics governing any 
given strange quark unaffected to a first approximation. In this case, one would expect 
the enhancement factor 1.115 to equally apply to the axial matrix element As. In view 
of the rough nature of this argument, a systematic uncertainty of 12% will be associated 
with this enhancement, thus covering the range of no enhancement of As up to twice 
the enhancement seen in the case of the scalar matrix element. 

Summarizing the diverse uncertainties and adjustments discussed above, the final 
estimates for the physical values of Jt s and As follows: 

/r. = 0.046(9)(l)(3)(5)(l) (35) 

where the uncertainties are, in the order written, statistical, due to finite lattice spacing, 
due to truncation of the chiral perturbation series, due to inadequate lattice dimensions, 
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and due to the adjustment of the strange quark mass to the physical value. In turn, 

A S = -0.031(16)(i?)(l)(4)(3)(4) (36) 



where the uncertainties are, in the order written, statistical, due to renormalization, 
due to finite lattice spacing, due to truncation of the chiral perturbation series, due to 
inadequate lattice dimensions, and due to the adjustment of the strange quark mass 
to the physical value. To quote a succinct final result, if one combines all the system- 
atic uncertainties discussed in this section together with the statistical uncertainty in 
quadrature, 



7 Conclusions 

This investigation focused on two of the most basic signatures of strange quark degrees 
of freedom in the nucleon, namely, the strange quark contribution to the nucleon mass, 
characterized by Jt s , and the portion of the nucleon spin contained in strange quark 
spin, As. A high amount of averaging not only of the disconnected strange quark loop, 
but also of the nucleon two-point function led to clear signals for both fx a and As 
especially at the heaviest pion mass, = 495 MeV, with the signals deteriorating, but 
not disappearing, as the pion mass is lowered to = 356 MeV and = 293 MeV. 
Combining all the lattice data, the signals survive chiral extrapolation to the physical 
pion mass. Systematic uncertainties remain under adequate control; the only source of 
systematic uncertainty which was not quantified is gluonic operator admixtures to fr a 
under renormalization. However, as discussed in section I4.1[ a scenario in which these 
admixtures rise to the level at which they begin to appreciably influence the conclusions 
reached regarding fx a seems highly implausible. The gluonic admixtures would have 
to be an order of magnitude larger than the related light quark admixtures, which 
were constrained to the 1% level. Nevertheless, a more quantitative corroboration of 
this argument would be desirable. All other systematic uncertainties were quantified, 
cf. section [6j and, while some of them are still sizeable, none rise to the level of the 
statistical uncertainties. With respect to controlling systematic uncertainties, the use of 
a (to a very good approximation) chirally symmetric discretization of the strange quark 
fields in the matrix elements (CO,© proved very advantageous, since it provides for 
benign renormalization properties, including the almost complete suppression of light 
quark admixtures to fx a alluded to above. This stands in contrast to, e.g., the case 
of Wilson fermions, in which the evaluation of fx a is considerably complicated by the 
presence of strong additive mass renormalizations [T2"||2"4]. 

The magnitudes obtained for both fx a and As appear natural. Neither quantity is 
abnormally enhanced; strange quarks contribute about 4.5% of the nucleon mass, and 



As 



0.046(11) 
-0.031(17) 



(37) 
(38) 
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the magnitude of the strange quark spin, which is polarized opposite to the nucleon spin, 
amounts to about 3% of the latter. The conditions provided by nucleon structure for 
dark matter detection via coupling of the Higgs field specifically to the strange quark 
component thus do not appear to be as favorable as assumed in the most optimistic 
scenaria. There is also no indication from the result for As that an unnaturally large 
contribution to the spin of the nucleon is hidden in the small-a; strange quark sector that 
has hitherto eluded experimental study. The strange quark spin does indeed appear to 
be polarized slightly in the direction opposite to the nucleon spin, as also indicated by 
the preponderance of phenomenological studies. However, the magnitude of As found 
in the present calculation is smaller than the magnitudes extracted in the analyses J5j|9] 
which assume a substantial enhancement of the strange quark helicity distribution at 
small momentum fraction x. 
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Appendix: Estimate of dm res j/dm h s are 

The quantity dm res ^/ dm h s &re enters the estimate of operator mixing effects in fx s , cf. sec- 
tion 14.11 No direct data for this quantity are available, but an order of magnitude es- 
timate can be constructed from related data on the residual mass obtained within the 
LHPC program and in the present work, summarized in Table |5j It should be noted 



amf are 


0.0081 


0.0313 


0.081 


10 3 • am reS) i 


1.6 


1.2 


0.7 


10 4 • am reS)S 


9.0 


8.1 


7.1 



Table 5: Residual masses m res i and m lress at varying tti 



bare 



that all these results were obtained at a constant lattice spacing a. Fitting parabolae to 
the data in Table [5] yields the following derivatives at the S77(3)-flavor symmetric point 
am s = ami = 0.081 and at the lightest mi considered in this work, ami = 0.0081: 



dm 



res, l 



-0.005 



dm 



res A 



am s =arni =0.081 



-0.020 (39) 



am 3 =0.081, am; =0.0081 
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dm r 



dm\ 



bare 



-0.0007 



dm rp 



am a =an%i =0.081 



5m bare 



-0.0045 (40) 



am 3 =0.081, ami =0.0081 



The simplest estimate for dm res j/dmg e from this can be obtained by noting that, at 
the S77(3)-fiavor symmetric point, 



dm 



res J 



am 3 =airii =0.081 



dm r 



9m bare 



-0.00035 



(41) 



=am i =0.081 



(where the factor 1/2 stems from the fact that d/drrii = d/dm u + d/drrid). Assuming 
that the derivative of m res> i in the m s -direction varies only weakly as one changes mi, 
one arrives at the estimate that also 



dm. 



res,l 



dm haie 



-0.00035 



(42) 



am 3 =0.081 , am; =0.0081 



Note that this is most likely an upper bound in magnitude, since one would expect the 
characteristics of the gauge fields entering the Dirac operator to be dominated by the 
light quark degrees of freedom relative to the strange quark degrees of freedom as mi 
becomes smaller. 

A check on this estimate can be constructed by the following alternative chain of rea- 
soning. First, note that the quantity dm res ^ldm\ we contains two contributions, namely, 
one from the explicit variation of mi in the Dirac operator which m reSi i characterizes, and 
the other from the implicit dependence of the gauge field ensemble on mi. To estimate 
dm res j/dm^ aTe , one therefore needs to apply two correction factors to dm res ^/ dmf^: 
A factor characterizing the proportion of the variation of m res j due specifically to the 
implicit variation of the gauge fields, and a factor characterizing the strength of that 
variation with m s as opposed to m/. The order of magnitude of these correction factors 
can be inferred as follows. The former factor is available at the 5'[/(3)-flavor symmetric 
point as the ratio 



dm reStS /dm\ 



bare 



dm reS) i/dm\ 



bare 



0.14 



(43) 



am a =arni =0.081 



For the purposes of the present argument, it will be assumed that this factor only varies 
mildly as mf are is lowered to amj 3are = 0.0081. On the other hand, assume also that 



dm r 



Q m b^e 



dm r 



implicit, ami =0.081 



dm°™ e 



(44) 



implicit, ami =0.0081 



i.e., the implicit variation of m reSyS via the dependence of the gauge fields on m s changes 
only mildly as a function of m^ Note that, while this assumption is analogous to the 
one leading to (14"2"|) . it is better founded since m reSjS itself varies less with mi than m Tes ^. 
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Again, one would expect the left-hand side to represent an upper bound for the right- 
hand side. Noting that the left-hand side is identical to one-half the quantities in the 
left-hand identity in (HOI) , one thus has 



q gg 0.0007/2 >.""(•<.«..-/ <."".s | linnlMi .null u.u^j (45) 

nnri/ic a / bare I v ' 



-0.0045 



<9m reSjS /<9m!? are 


implicit, ami =0.081 


<9m reSjS /<9mj? are 


am ; =0.0081 



dm res , s /dm 



bare I 



| implicit 



dm res ^ s /dmf arc 
dm res ,i/dm 



bare 



dm reS: i/dm h£L " 



l | implicit 



ami =0.0081 



am;=0.0081 



(46) 
(47) 



where in the final step it is assumed that changes in the numerator and denominator 
due to changing the quark mass in the Dirac operator from m s to mi approximately 
cancefl. This is the desired conversion factor characterizing the strength of the implicit 
variation of m res j with m s relative to the one with mi. Correcting, as proposed above, 
dm res j/dm]' a ' re by the two factors (|4*3|) and ( |4"T|) . one finally arrives at the alternative 
estimate 



dm. 



resl 



« 0.08 • 0.14 • (-0.020) = -0.0002 , (48) 

am a =0.081,am; =0.0081 



consistent in order of magnitude with (142]) . especially in view of the latter being expected 
to represent an overestimate. 

The estimate (|4*2l) for dm reS) i/dm^ are at the lowest light quark mass considered in 
the numerical calculations in this work is used in section 14.11 to constrain the influence 
of operator mixing effects in the renormalization of fx s - 
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